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We point out that the potential of Goldberger and Wise for stabilizing the distance between two 
3-branes, separated from each other along an extra dimension with a warp factor, has a metastable 
minimum when the branes are infinitely separated. The classical evolution of the radion (brane 
separation) will place it in this false minimum for generic initial conditions. In particular, infla- 
tion could do this if the expansion rate is sufficiently large. We present a simplified version of 
the Goldberger- Wise mechanism in which the radion potential can be computed exactly, and we 
calculate the rate of thermal transitions to the true minimum, showing that model parameters can 
be chosen to ensure that the universe reaches the desired final state. Finiteness of bulk scalar field 
brane potentials can have an important impact on the nucleation rate, and it can also significantly 
increase the predicted mass of the radion. 
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I. INTRODUCTION 

One of the most striking proposals in current elemen- 
tary particle theory is the existence of extra dimensions 
which arc hidden from us, not necessarily by their small- 
ness, but by our confinement to a four-dimensional slice 
(a 3-brane) of the full spacetime [1]. Randall and Sun- 
drum [2] have produced a verion of this scenario which 
is particularly attractive because of its apparent links 
to deep theoretical developments: the conformal field 
theory/5D-anti-de Sitter correspondence and holography 
[3]. On a more practical level, their idea provides an ex- 
planation for why the Planck scale is so much greater 
than the weak scale, independently of supersymmetry. 
The cosmological implications of this model have also 
been the subject of vigorous study [4]. 

The Randall- Sundrum proposal involves a "Planck 
brane" located at a position y = in a single additional 
dimension, and a second "TeV brane" located at y = 1, 
in our conventions. The extra dimension is permeated by 
a negative bulk vacuum energy density, so that the space 
between the branes is a slice of anti-de Sitter space. Solv- 
ing the 5-D Einstein equations results in the line element 

ds 2 = e -2<r(y)( dt 2 _ dS ^ _ b 2 dy 2. ^ 

a(y) = kby, y S [0, 1] (2) 

The constant k is related to the 4-D and 5-D Planck 
masses, M and M p respectively, by 



where M p 2 = 8nG N and the 5-D gravitational action is 
S = |M 3 / d 5 Xy/^R. The warp factor e~ CT ^ deter- 
mines the physical masses of particles on the TeV brane: 
even if a bare mass parameter mo in the TeV brane La- 
grangian is of order M p , the physical mass gets rescaled 
by 
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Wphys = er a m n ; a = a(l) = kb, (4) 

as can be easily derived starting from a scalar field 
action written covariantly in terms of the metric (1), 
S = | J d 4 xe~ 4CT (e 2<T ((90) 2 — m§0 2 ), and rescaling 4> so 
that the kinetic term becomes canonically normalized. 

In order for m to be of order 100 GeV, the combina- 
tion kb must be approximately 36, so that e~ CT ~ 10 -16 . 
Yet in the original proposal, the value of b, which is 
the size of the extra dimension, was completely undeter- 
mined. It is a modulus with no potential, which is phe- 
nomenologically unacceptable. For one thing, the par- 
ticle associated with 4-D fluctuations of b, the radion, 
would couple to matter on the TeV brane similarly to 
gravitons, but more strongly by a factor of e CT [5]. This 
would lead to a fifth force which could easily have been 
detected. Furthermore, a massless radion leads to prob- 
lems with cosmology: our brane-universe would have to 
have a negative energy density to expand at the expected 
rate, assuming that energy densities on the branes are 
tuned to give a static extra dimension [6] . It was shown 
in refs. [7,8] (see also [9]) that this problem disappears 
when the size of the extra dimension is stabilized. 

Radion stabilization is therefore a crucial ingredient of 
the Randall-Sundrum idea. Goldberger and Wise have 
presented an elegant mechanism for accomplishing this 
[10], using a bulk scalar field. Self-interactions of the 
field on the branes forces it to take nonvanishing VEV's, 
Vq and V\ respectively, which are generally different from 
each other. The field thus has a gradient in the extra 
dimension, and the competition between the gradient and 
potential energies gives a preferred value for the size of 
the extra dimension. In other words, a potential for the 
radion is generated, which has a nontrivial minimum. It 
is easy to get the correct brane separation using natural 
values of the parameters in the model. 

Roughly speaking, the radion potential has the form 
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v{<f>) - \<p 4 



Vl 



(5) 



with a nontrivial minimum at <j> = /(wi/w ) 1/e - How- 
ever, there is another minimum at <fi = 0, which repre- 
sents an infinitely large extra dimension. This could not 
describe our universe, because then e _CT would be zero, 
corresponding to vanishing particle masses on the TeV 
brane. In the more exact expression for the potential, we 
will show that <p — is actually a false vacuum — it has 
higher energy than the nontrivial minimum. Neverthe- 
less, is is quite likely that the metastable state could be 
reached through classical evolution in the early universe. 
The question then naturally arises whether tunneling or 
thermal transitions to the desired state occurs. This is 
the subject of our study. 

Such detailed questions about the viability of the 
Goldberger-Wise mechanism are important because there 
are few attractive alternatives at the moment. Rcf. [11] 
recently studied Casimir energies of fields between the 
branes as a possible origin of a stabilizing potential. They 
found that stabilization in this way is indeed possible, 
but that the resulting radion mass is too small to be 
phenomenologically consistent if the size of the extra di- 
mension is that taken to be that dictated by the hierarchy 
problem. 

In section 2 we derive the Goldberger-Wise potential 
for the radion in a slightly simplified model which al- 
lows for the exact solution of the potential. The classical 
evolution of the radion field is considered in section 3, 
where we show that for generic initial conditions, the 
universe reaches a state in which the radion is not stabi- 
lized, but instead the extra dimension is expanding. This 
is a metastable state however, and in section 4 the rate 
of transitions to the minimum energy state, with finite 
brane separation, is computed. Conclusions are given in 
section 5. 



II. RADION POTENTIAL 

Let i/j(y) be the bulk scalar field which will be respon- 
sible for stabilizing the radion. Its action consists of a 
bulk term plus interactions confined to the two branes, 
located at coordinate positions y = and y = 1, respec- 
tively. Using the variable a of eq. (2) in place of y, the 
4-D effective Lagrangian for a static ip configuration can 
be written as 
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((d a iP) 2 + mV) da 



woOo - v ) - e 4CT rai(V>i - vi) 



(6) 



where rh is the dimcnsionless mass m/k, ipi are the values 
of ip at the respective branes, and the orbifold symmetry 
ip(a) — ip(—(j) is to be understood. In 5-D, the field ip, as 



well as the VEV's Vi on the two branes, have dimensions 
of (mass) 3 / 2 , while the parameters rrii have dimensions 
of mass. Denoting d a ip — ip' , the the Euler-Lagrangc 
equation for ip is 

ke~ ia {ip" - Alp' - m 2 iP) = 2m (ip - v )S(cr) 

+ 2e- 49 m 1 (ip 1 - Vl )S(a-a). (7) 

The general solution has the form 

ip = e 2a (Ae va + Be-") ; 



2 + e. 



(8) 



To get the correct hierarchy between the Planck and weak 
scales, it is necessary to take rh 2 to be small, hence the 
notation e. 

The brane terms induce boundary conditions specify- 
ing the discontinuity in ip' at the two branes. Imposing 
the orbifold symmetry ip(—<r) = ip(a), this implies that 



^'(0) = rh (ip - v ) 
ip'(a) = -rai(V>i - vi), 



(9) 



where we defined rhi — nii/k. Substituting the solution 
(8) into (9) allows us to solve for the unknown coefficients 
A and B exactly. In this respect, the present model is 
simpler than that originally given in ref. [10]. There the 
brane potentials were taken to be quartic functions, so 
that A and B could only be found in the approximation 
where the field values ipi were very strongly pinned to 
their minimum energy values, vi. In our model this would 
occur in the limit m, — > oo. However, we can easily 
explore the effect of leaving these parameters finite since 
A and B can be determined exactly. Let us denote 



= e 



-kb 



= e 



(10) 



which will be convenient in the following, because <p is 
proportional to the canonically normalized radion field. 
Then A and B are given by 
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-d<p v + C 2( p 2 )<p v /D(<P) 



b= (< ', " -r. 



2 ) 4> V /D(4>) 
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where 



D{<t>) 



Ci = m v (mi - e) 
C 2 = mivi(m + e) 
C 3 = m v (rhi + 4 + e) 
C4 = rhiVi (rriQ — 4 — e) 

(C 2 C 3 - CiC 4 2 » 



(m miv vi) 



(12) 



It can be checked that in the limit rhi — ► oo, the field 
values on the branes approach ipi — > Vi. For finite rhi, 
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the competing effect of minimizing the bulk energy causes 
departures from these values, however. 

The solution for can be substituted back into the 
Lagrangian (6) to obtain the effective 4-D potential for 
the radion, V{4>). However, rather than substituting di- 
rectly, one can take advantage of the fact that ip is a 
solution to the Euler-Lagrange equation. Doing a partial 
integration and using the boundary terms in (7) gives a 
simpler expression for V(cj>). In the general case where 
the brane potentials are denoted by Vi(ip), one can easily 
show that 



(13) 



In the present case, we obtain 

v{4>) = -c 

= mov Q (vo - ip ) + ^miv^vi - V>i) 
= m v (v - (A + B)) 

+ $ 4 m 1 v 1 (v 1 -4>- 2 (A4>- v + 4>"B)) (14) 

In the following, we will be interested in the situation 
where V(<f>) has a nontrivial minimum for very small val- 
ues of <j) ~ 10~ 16 , as needed to address the weak scale 
hierarchy problem. It is therefore a good approximation 
to expand V(<f>) near <p = 0, keeping only the terms which 
are larger than 4> 2lJ . This is accomplished by expanding 
the denominator D{<p) in eqs. (11), after which the po- 
tential can be written in the simple form 

V(4>) = A0 4 ((1 + e 4 - £l ) ^ - 2r, (1 + e 4 ) 4> e + ry 2 ) 

(15) 

where we introduce the notation 

Q = 7; £o = ei = (16) 

4 mo mi 



ri= (1 + co) 



Vl 

v 



and 



( 1 + ^T+ £ o) 



(17) 



(18) 



In the following it will be convenient for us to rewrite 
V((f>) in the form 



where c± are given by 

'(l + e 4 )± V /(l + e 4 ) 2 -(l + e4-£i) 



C± = 7] 



(1 + £4 - ei) 



(19) 



(20) 



and A' = A(l + e 4 - ei). 



III. PHENOMENOLOGY AND EARLY 
COSMOLOGY OF THE MODEL 

The warp factor which determines the hierarchy be- 
tween the weak and Planck scales can be found by mini- 
mizing the potential (19). Expanding in e, it has a global 
minimum and a local maximum at the respective values 
<t>+ and </>_: 



The last approximation holds in the limit of small e, e 
and ei; it is not always an accurate approximation for 
the parameter values of interest, so we will use the exact 
expression in any computations which might be sensitive 
to the actual value. The positions of the zeroes of V, 
<f) = C± e , are slightly greater than <p±, by the factor (1 + 



e 4 ) 



iA 



=1/4 



as can be seen by comparing (21) with 



(20). The large hierarchy of <fi+ ~ 10 -16 is achieved by 
taking a moderately small ratio v\/vq < 1 and raising 
it to a large power,* 1/e. This leads to the mass scale 
which functions like the cutoff on the TcV brane, 



-Mr, = M' 



TcV 



(22) 



where Mx e v/(1 TeV) is a number of order unity, which 
we will take to be one of the phenomenological free pa- 
rameters of the model. The choice of Mt c v specifies pre- 
cisely where we want our cutoff scale to be. In ref. [10] 
the exponential corrections to (vi/vq) 1 ^ in (21) were not 
considered; these change somewhat the values one might 
choose for Vi/vq and e to get the correct hierarchy. The 
factor e ± v / ( 1 / mi+1 / 4 )/ e in particular can be significant. 

Refs. [7] and [5] showed that the canonically normal- 
ized radion field is <j> — f<j>, where / = y 6M 3 /fc is an- 
other scale of order M p J The 4-D effective action for the 
radion and gravity is 
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: (y%4>d»$-v{4>j) 



(23) 



*An alternative possibility, taking vi/vo > 1 and e < 0, cor- 
responding to a negative squared mass in the bulk Lagrangian 
(6), does not work. Although the negative m 2 does not neces- 
sarily lead to any instability, since the field is prevented from 
going to infinity by the potentials on the branes, the radion 
potential has no non trivial mi nimum in this case. 

^The choice / = y // 24M' i /k in ref. [5] seems to correspond 
to an unconventional normalization for M„. 
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where V(<j>) is the potential (19) and R is the Ricci scalar. 
To get the correct strength of gravity, we must therefore 
have 



M 3 



M, 
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f = V6M P . 



(24) 



The radion mass is / 2 times the second derivative of V 
evaluated at its minimum. We find the value 

ml = 4 A (l + | - ^ (VTT^ + 2j£) V 2 ft e 3 / 2 , 

(25) 

which implies that m^ is of typically of order e 3 / 4 times 
the TeV scale. The factor of e 3 / 4 leads to the prediction 
that the radion will be lighter than the Kaluza-Klein ex- 
citations of the graviton, which would also be a signal 
of the new brane physics [5]. However, we see that the 
corrections due to finite mi, which were not explicitly 
considered in [5] , can possibly compensate this and make 
the radion heavier, if rh\ ~ e. 

Now let us turn to the evolution of (f> in the early uni- 
verse. For this purpose it is important to understand 
that the depth of the potential at its global minimum, as 
well as the height of the bump separating the minimum 
from (j> = 0, is set by the TeV scale. The values of the 
potential at these extrema are approximately (to leading 
order in 64, but exact in ei) 



£4 



1 + £4 - ei 



(l±v^)- 



(26) 



Since A - M 4 



the depth at the minimum is V{4>+) ~ 
-e 3 / 2 0(TeV) 4 — suppressed slightly by the factor of e 3 / 2 . 
The height of the bump at <fi- can be considerably smaller 
because of the exponential factors in (21). In fact 



V{<t>. 



where 



1 — \/ei + e4 




(27) 



exp(-^a + Sr)) (28) 



For example, if e = 0.01 as suggested by [10], ft 4 is 
less than 10~ 17 . If the brane potential parameter mi is 
not large, so that mi £ 1, the suppression will be much 
greater. Figure 1 illustrates the flatness of the potential 
for the case e = 0.2, where the barrier is not so sup- 
pressed. The new mass scale f2 TeV <C 1 TeV is due to 
the small curvature of the radion potential at the top of 
the barrier, and its smallness will play an important role 
in the following. 

Thus the barrier separating the true minimum at 
from the false one at <j> — is extremely shallow. More- 
over a generic initial condition for the radion is a value 



like </> ~ 1, quite different from the one we want to end 
up with, <f> <~ 10 -16 . Clearly, the shape of the potential 
is such that, if we started with a generic initial value for 
4> 1 it would easily roll past the local minimum and the 
barrier, hardly noticing their presence. The point = 
toward which it rolls is the limit of infinite brane sep- 
aration, phenomcnologically disastrous since gravity no 
longer couples at all to the visible brane in this limit. 

One might wonder whether inflation could prevent this 
unwanted outcome, since then there would be a damping 
term in the (j) equation of motion, possibly causing it to 
roll slowly: 



$ + 3H4> + V^(<f>)=0. 



(29) 



Indeed, with sufficiently large Hubble rate iJ, the mo- 
tion could be damped so that <j> would roll to its global 
minimum. The condition for slow-roll is that 



Kff « 9ff 2 



(30) 



However, inflation is a two-edged sword in this instance, 
because the effective potential V e s for the radion gets ad- 
ditional contributions from the curvature of the universe 
during inflation. From eq. (23) one can see that 



v eS (4>) = v($) 



M 3 

~2k ^ 

M 3 



= V{<j>)+H 



k 

2i2 



(31) 



using the relation R — 12H 2 which applies for de Sitter 
space, and eq. (24). The new term tends to destroy the 
nontrivial minimum of the radion potential. One can 
estimate the relative shift in the position of the minimum 

as 



s$+ sv'(4> + ) _i? 2 

" ~0 + V"(0+) " m% 



(32) 



This should be less than unity to avoid the disappearance 
of the minimum altogether. 

Combining the requirement that the global minimum 
survives with the slow-roll condition (30), evaluated near 
the minimum of the potential, we find the following con- 
straint on the Hubble rate: 



\m\ < H 2 < ml. 



(33) 



This is a narrow range, if it exists at all. In fact, one 
never expects such a large Hubble rate in the Randall- 
Sundrum scenario since the TeV scale is the cutoff: H 
should never exceed T 2 /M p - 10~ 16 TeV if the classical 
equations are to be valid. The problem of radion stabi- 
lization might also be exacerbated because contributions 
to the energy density of the universe which cause infla- 
tion can give additional terms of the type (j) 2 to V e s which 
are not considered in the above argument. For example, 
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Figure 1: V(<f>) versus <f> for e = 0.2, showing the smallness of the barrier (right) relative to the minimum (left). Notice the 
difference in vertical scales. 



a field in the bulk which does not have its endpoints fixed 
on the branes has a 5-D energy density which is peaked 
on the visible brane [12], ps{y) ~ poe 2kby , and gives a 
contribution to V c g of 

6V ~ po f dybe- ikhy p e 2kby 
Jo 

= bp 4> 2 , (34) 

remembering that e~ kb = <j>. Such a contribution could 
destroy the nontrivial minimum even if (33) is satisfied. 
In any case, it does not appear to be natural to tune the 
Hubble rate during inflation to try to solve the stabiliza- 
tion problem. 



and consider the transitions of the radion starting from 
the time of reheating. The criterion that the phase tran- 
sition completes is that the rate of bubble nucleation per 
unit volume, T/V, exceeds the rate of expansion of the 
universe per Hubble volume: 

£ £ H A (35) 

The reason is that the bubbles expand at nearly the speed 
of light, so the relevant volume is determined by the 
distance which light will have traveled by a given time, 
which is of order l/H. 

A. The Euclidean Bounce 



IV. PHASE TRANSITION TO THE TRUE 
VACUUM 

Since the barrier of the radion potential is too small 
to prevent the radion from rolling into the false mini- 
mum, perhaps we can take advantage of this smallness 
to get tunneling or thermal transitions back into the true 
vacuum. The situation is quite similar to that of "old in- 
flation" [13], except that in the latter, the transition was 
never able to complete because the universe expanded too 
rapidly compared to the rate of nucleation of bubbles of 
the true vacuum. In the present case this problem can 
be avoided because we are not trying to use the radion 
for inflation. Indeed, a small amount of inflation may 
take place before the tunneling occurs, since the radion 
potential is greater than zero at (f> = 0, but we will not in- 
sist that this be sufficient to solve the cosmological prob- 
lems inflation is intended to solve — otherwise we would 
be stuck with the problems of old inflation. Instead we 
will assume that inflation is driven by some other field, 



To compute the nucleation rate T/V, one must con- 
struct the bounce solution which is a saddle point of the 
Euclideanized action [14], in other words, with the sign of 
the potential reversed. This is a critical bubble solution 
with the boundary conditions 

<p(r)\ r =o = <Po; (p'(r)\ r =o = 0; 
4>{r)\ r ^oo -0; 0'(r)|™ = 0. (36) 

The value of (f> n which ensures the desired behavior as 
r — > oo cannot be computed analytically because the mo- 
tion of the field is damped by the term <j>'/r in the equa- 
tion of motion. We will consider bubble nucleation at 
finite temperature in the high T limit, where the bounce 
solutions are three dimensional. The equation of motion 
is 

1 (rV)' = +VM, (37) 

where now V e s includes thermal corrections, which are 
much larger than the H 2 (j) 2 term considered in eq. (31): 
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Figure 2: The bounce solution, for the parameters e : 
4> — > <j> and r — > f is given in eqs. (43-44). 



0.01, mi = 0.1, m 4 , = T = 100 GeV and M TcV = 1 TeV. The rescaling 



Cb 



V eS {<t>) 

m 2 h (0) = V^) + ^T 2 ^ 4 >(O) 



(38) 
(39) 



where q, = 3.9 if the renormalization scale is taken to 
be equal to T. The leading thermal correction is of or- 
der T 2 (j) 2 , whereas H 2 4> 2 is suppressed by T 2 /M 2 rela- 
tive to this. The cubic term m\ h becomes imaginary if 



V"(<j>) 



24 



V"( 4 )(0) becomes negative; we take the real 



part only. The term j^T 2 V^ (0) represents the thermal 
mass, which appears in the cubic and quartic terms when 
ring diagrams are resummed [15]. We subtract a constant 
term V from V e fi so that V^fj(O) = 0, as is needed to 
properly compute the action associated with the bounce 
solution. 

The thermal corrections to the effective radion poten- 
tial cause the bounce solution to fall exponentially at 
large r: 



£ p-r/ro 



(40) 



where l/r ~ r)y/A/f 4 T if A// 4 < 1 [the exact ex- 
pression is 1/ro = \fJjJ*T, in terms of quantities to be 
defined below, in eqs. (43) and (55)]. Once the bounce 
solution is known, it must be substituted back into the 
action, which can be written as 



4^ 

Y 



drr 2 (^ /2 + V cS ^)). 



The nucleation rate is given by 
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s) '^" 1/2e " s 



(41) 



(42) 



where w_ is the imaginary frequency of the unstable 
mode of fluctuations around the the bounce solution, and 



V is the fluctuation determinant factor, to be described 
below. A typical profile for the bounce solution is shown 
in figure 2. 

For the numerical determination of the bounce solution 
and action, as well as understanding their parametric 
dependences, it is convenient to rescale the radius and 
the field by 



VAT' 
ZT~4>- 



A- /4 c_ 



fc 



Then the action takes the form 



S=—S(e,m 1 ,\,Z); 



(43) 
(44) 

(45) 
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J^drf 2 ^^ 2 + 4> 2 f 2 ) 



+ z 2 4> 4 f 



i Va 

12tt Z 2 



+ 12Z 



where 



fS) = & - m e - 

k(« = (i + ^i)/»W 



3/2' 



(46) 



(47) 



In the following, it will be helpful to keep in mind that 
Z can be extremely small, of order Vi in (28) when e is 
small, whereas V\ tends to be closer to unity, depending 
on the mass of the radion and the definition of the TeV 
scale (22): 
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2^6 M Te v 



(1 — v^T+ei) 
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4 + 2^") 
mi y 



1/2 



(48) 



1 
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In the rescaled variables, the value (f> = 1 corresponds 
to the first zero of the potential, which would be the 
starting point of the bounce if energy was conserved in 
the mechanical analog problem, i.e., if there was no vis- 
cous damping term <p' /r in the equation of motion. The 
actual starting point turns out to have a value in the 
range 4>q ~ 1.5 — 3 because of this. The rescaled action S 
depends mainly on the model parameters e and mi, for 
the parameter values which are of interest to us. All the 
sensitive exponential dependence on e, namely the factor 
c^ e , is removed from S. Numerically we find that 



fmi(emi)- 3 / 4 , 



(49) 



except when rh\ becomes close to e. For mi slightly 
smaller than e, ei starts to exceed 1, and c_ becomes 
negative, signaling the onset of an instability in the ra- 
dion potential toward coincidence of the two brane po- 
sitions. Figure 3 shows the dependence of log(,S/mi) 
versus log(emi). 
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Figure 3: \og w (S/rh\) versus log 10 (emi), where S is the 
rescaled bounce action, eq. (46). The other parameters are 
m = T = 100 GeV and M To v = 1 TeV. 

We have computed the bounce solution and the cor- 
responding action for a range of parameters e and rh\ 
which can be consistent with the solution to the hierarchy 
problem (i.e., that <j> <~ 10~ 16 at the global minimum). 
The size of the bounce in position space, measured as 
the width at half- maximum, is small near e = mi, and 
reaches a larger constant value as mi — ► oo. Using the 
rescaled radial variable f — r^/XT, the dependence of the 
width on e and rh\ is shown in figure 4. 
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Figure 4: Half-width w of the bounce solution, in terms of 
the rescaled radial distance f — rvAT, versus log 10 (rai), for 
the same parameters as in figure 3. 

The action of the bounce can be much greater than or 
much less than 1, depending on the parameters: for e ~ 
mi < 1, S < 1, while for larger values of e and mi, S » 
1. Where the crossover occurs (S <~ 1) depends on m^, 
T and MxeV- This behavior can be inferred from figure 3 
(showing S) and the dependences of the coefficient in the 
relation S = (Z 2 /\/X)S. Rather than presenting further 
results for S directly however, we will turn to the more 
relevant quantity, the rate of bubble nucleation. For this 
we need to determine the prefactor of e~ s in the rate. 



B. Prefactor of Bubble Nucleation Rate 

The bounce action is the most important quantity de- 
termining the rate of tunneling, since it appears in the 
exponent of the rate (42). Since we do not have a model 
for the inflation and reheating of the universe which must 
occur prior to the bubble nucleation, hence an exact pre- 
diction for the reheating temperature which enters the 
rate, it would not be worthwhile to compute the prefac- 
tors in eq. (42) very accurately; however we can estimate 
their size. 

First, consider the frequency of the unstable mode. 
uP_ is the negative eigenvalue of the Schrodingcr-likc 
equation for small fluctuations 8<j> around the bounce so- 
lution, which we will denote by <j>b(r): 



+ 



+ 



d 2 V cS 



oj 2 _5(j) (50) 



Mr) 



Rescaling the radius and background field exactly as in 
eqs. (43-44), eq. (50) becomes 



(51) 
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Figure 5: The potential for small fluctuations around the bounce solution, U(f) = X 2 V" (4>b(f)), as a function of f, for the 
parameters e = 0.02, mi = 5 (left) and e = 0.1, mi = 25 (right). 



[/(f) 



1 d 2 V eS/7 



AT 2 



+ ^A(^/ 4 + X(/ 4 + 2/|// 2 )) 



12(^/a) 2 



7& 



(52) 



where primes now denote J=, X = 12Z 2 2 / 2 , and the fi 
are defined in eq. (47). Except when the radion mass is 
signficantly larger than 100 GeV, A is much smaller than 
1, and U is dominated by the first two terms in (51). Of 
these, the first (/ 4 ) always dominates if Z <C 1, while 
the second (X) can be important near r = if Z > 1. 
The two different cases are illustrated in figure 5. Since 



1 < ( 
that 



< — , both fi and X are negative at r = 0, so 



U = U(0) £* -e 3 /Vl + £ 



(2 + ln^o + 12Z 2 <^(i+ln 



0o)) (53) 



and thus the smallest eigenvalue of eq. (50) is negative. 
This is the unstable mode of the saddle point solution, 
with imaginary frequency of order 

u?_ ~ U XT 2 . (54) 

Recall that \ui- | appears in the prefactor of the nucleation 
rate T/V. 

As f —* oo, U(f) approaches a maximum value 
U. = U(oo) = £ (l - f ^ + ) , (55) 

which determines the asymptotic behavior of the fluctu- 
ations at large f: 5<p ~ e~ v/ ^ 7 " r . The fluctuations around 
the false vacuum state (</> = 0) thus have a mass given 
by 



m 



(56) 



which will be relevant for the following. 

Next we must estimate the functional determinant fac- 
tor, 



V 



det'(-d 2 - V 2 + V'jfo)) 
det(-<9 2 - V 2 + V"(0)) 



(57) 



where r is imaginary time (t e [0, 1/T]), V 2 is the three- 
dimensional Laplacian, </>b is the bounce solution, and 
the prime on det' means that the three translational 
zero-mode eigenvalues must be omitted from the deter- 
minant for fluctuations around the bounce. These zero 
modes correspond to spatial translations of the bubble 
solution. Because of their removal, T> has dimensions of 
(energy) -6 , as is required to get a rate per unit volume 
in eq. (42). 

Ref. [16] has given a thorough account of how to 
compute I? by a method which was discussed for one- 
dimensional systems in [17]. In 3-D one should classify 
the eigenvalues of the fluctuation operators by the quan- 
tum numbers n and I, denoting Matsubara and angu- 
lar momentum excitations, respectively. Then T> can be 
written as a product, V — \\ n l V n j. 

Ref. [16] shows that the contribution to V from the 
Zth partial wave can be expressed, to leading order in a 
perturbative expansion in the potential U(r), as 



21+1 



(i+h^y . (58) 

The quantity has the Green's function solution 

f-OC 

=2 / drrI l+1/2 {Kr)K l (Kr) {V {^{r)) - m 2 ) 
Jo 

t>OC 

= 2/ dffl^i^Kd^) (U(f)-U.) (59) 

J 
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using the modified Bessel functions / and K, and the 
mass m of the field in the false vacuum, eq. (56). For 
general Matsubara frequencies, v — 2imT, one defines 
k = \Jm 2 + v 2 . 

The subdeterminant for the n = (zero-temperature) 
sector of the theory has the usual ultraviolet divergences 
of quantum field theory, namely the vacuum diagram (y 
(the dot represents one insertion of V(<p)), which should 
be absorbed by renormalization of the zero of energy for 
the radion potential. Since we are not attempting to solve 
the cosmological constant problem here, we are going to 
ignore all of this and compute only the factor 2\i> which 
contains the translational zero modes — or more precisely, 
which has the zero modes removed. This removal is ac- 
complished by replacing 



l + h\ 



(i) dh\ 



(i) 



dn 2 



(60) 



Notice that this quantity has dimensions of (mass) -2 , 
and there are 21 + 1 = 3 such factors, so that |2?|~ 1//2 has 
dimensions of (mass) 3 , as required. From eq. (59) one 
can show that 



dK 2 
Iu = 



■Iu, 



\T 2 U 2 

*dyy 2 (VMtfifo))' (uijfc) u.) . (61) 



We have numerically evaluated the integral I v for each 
set of parameters. Our estimate for the fluctuation deter- 
minant factor in the nuclcation rate can then be written 



as 



IXH-V2 



XT 2 U 2 \ 3/2 
Iv 



(62) 



C. Results for Nucleation Rate 

Putting the above ingredients together to find the rate 
of bubble nucleation per unit volume, Y/V, we see that 
the latter depends on five undetermined parameters: e, 
rhi, m$, MxeV and the temperature T. Ref. [7] showed 
that, as long as the energy density on the TeV brane 
is much less than M| cV , the usual 4-D effective theory 
governs the Hubble rate: 



H 2 



P 



(63) 



where p is the total energy density, 



71 rpA 

5 *30 

U 2 + v{4>) - v(<j>+) 



±M 2 rn- 



(64) 



We take the number of relativistic degrees of freedom, 
<7», to be 100. The kinetic energy of the radion is 
zero since <p — in the false vacuum, so p^ is essen- 
tially the potential energy of the radion in the false vac- 
uum, assuming the 4-D cosmological constant is zero: 
V{4>) - V{4>+) = \V(4>+)\, which is given by eq. (26). 
Depending on the parameters, this can be comparable in 
size or dominate over the energy density of radiation. Us- 
ing our estimates for the prefactor of the tunneling rate, 
the logarithm of the ratio ofT/V to H 4 can be written 
as 



In 



VH 4 



In 



XWUoU, 
(2tt) 5 / 2 



S 



3/2 



(65) 



where S is the action of the bounce solution. The cri- 
terion for completion of the phase transition to the true 
vacuum state is that \n(Y /VH 4 ) > 0. The saddle point 
approximation leading to eq. (42) is only valid if the ac- 
tion S is not much less than 1. Otherwise, the barrier 
is not effective for preventing the field from rolling to 
the true minimum, as in a second order phase transition. 
This situation occurs in the vicinity of ln(Y /VH 4 ) <~ 150 
in the following results; thus the transition region where 
Y /VH 4 = 1 is well within the realm of validity of the 
approximation. 

In figure 6 we show the contours of constant 
ln(r /VH 4 ) in the plane of log 10 (mi) and log 10 (e), start- 
ing with the fiducial values T = = 100 GeV, Mt c v = 
1 TeV for the other parameters, and showing how the re- 
sults change when any one of these is increased. The de- 
pendences can be understood from the prefactor Z 2 /\f\ 
in the action, eq. (45): 



~ e 



3 /4_M TeV fi 2 



rp2 



(66) 



where we recall that Z and A are given by (43-44) and fi 
by (28). The factor SI is responsible for suppressing the 
bounce action when t < 1 or erh\ <C 1, explaining the 
shape of the allowed regions in each graph. Nucleation 
of bubbles containing the true minimum becomes faster 
when the temperature or the radion mass is increased, 
but slower if the definition of the TeV scale in increased. 
These dependences are dictated not only by the size of 
the barrier between the two minima in the effective po- 
tential, but also by the size of the bubbles. 

Interestingly, the borderline between allowed and for- 
bidden regions of parameter space falls within the range 
which is relevant from the point of view of building a 
model of radion stabilization. That is, some choices 
which would otherwise have been natural and acceptable 
are ruled out by our considerations. We sec furthermore 
that the choice of mi — > oo, as was effectively focused on 
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Figure 6: Contours of ln(r/V.ff 4 ) in the plane of log 10 (rhi) versus log 10 (e). The shaded regions are where the tunneling rate is 
too small for the phase transition to complete. Figure (a) has T = = 100 GeV, MtcV = 1 TeV. The other figures are the 
same except for the following changes: (b) T = 400 GeV; (c) m = 400 GeV; (d) M TcV = 2 TeV. 



in ref. [10], is not the optimal one for achieving a large 
nucleation rate. 

It might be thought that our analysis is rendered less 
important by the fact that one can always obtain fast 
nucleation simply by going to high enough temperatures. 
However it must be remembered that the TeV scale func- 
tions as the high-energy cutoff in the Randall-Sundrum 
scenario: the whole semiclassical description breaks down 
at super- TeV scales, where quantum gravity effects start 
to become important. From this point of view, the tem- 
peratures of 100 — 300 GeV which we are discussing are 
already rather high, and a fairly efficient mechanism of 
reheating at the end of inflation will be needed to gener- 
ate them. 



V. DISCUSSION 

In this paper we have presented a somewhat simpler 
model of radion stabilization by a bulk field (tp) than 
that of Goldberger and Wise [10]; although the physics 
is qualitatively identical, we are able to write the radion 
potential exactly, and thus explore the effect of letting the 
stabilizing field's VEV's on the branes be pinned more or 
less strongly to their minimum energy values. One such 
effect is that the mass of the radion can be significantly 
increased for small values of the parameter mi, which is 
the coefficient of the potential for ip on the TeV brane. 
Moreover if m\/k = mi is accidentally close to e, ap- 
proximately the minimum value consistent with a stable 
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potential, the radion mass can start to diverge, by the 
factor (1 — 4 — ^—)~ 1 / 2 . This modifies somewhat the ex- 

V 4 mi / 

pectation expressed in rcf. [5] that the radion mass will 
be small relative to the TeV scale, due to a factor of e 3 / 4 . 

Our main focus was on the problem that the radion po- 
tential has a local minimum at infinite brane separation, 
and that the barrier between the true and false minima 
is so small that for generic initial conditions, one would 
expect the true minimum to be bypassed as the radion 
field rolls through it. We showed that for a large range of 
parameters, the high-temperature phase transition to the 
true minimum is able to complete, thus overcoming the 
problem. There are however significant constraints on 
the model parameters, and the initial temperature after 
inflation, to insure this successful outcome. 

There remain some outstanding issues. The form of 
the radion effective potential is such that the field is able 
to reach <f> = in a finite amount of time; yet = 
represents infinite brane separation in the extra dimen- 
sion. This paradoxical situation may be due to the as- 
sumption that the stabilizing field, Vj is always in its 
minimum energy configuration at any given moment. In 
reality tp must require a finite amount of time to respond 
to changes in the radion. Thus one should solve the cou- 
pled problem for time- varying <f> and tp to do better. This 
is probably a difficult problem, which we leave to future 
study. 

A related question is whether it is correct to treat ther- 
mal fluctuations of the radion field (f> analogously to a 
normal scalar field with values in the range (—00,00). 
Since <p is related to the size of the extra dimension by 
4> = fe~ kb , its range is [0, /]. We have not studied what 
effect this might have on the thermal part of the effective 
potential; instead we assumed that the usual treatment 
suffices. 

Another approximation we made was to ignore the 
back-reaction of the stabilizing field on the geometry 
Ref. [18] has given a method of finding exact solutions to 
the coupled equations for the warp factor a(y) and the 
stabilizing field ip(y). This method cannot be applied 
in the present case because it works only for bulk scalar 
potentials with a special form that, among other things, 
requires them to be unbounded from below. * Moreover, 
since the method of [18] generates only static solutions to 
the equation of motion, it cannot be used to deduce the 
radion potential, which is a probe of the response of the 
geometry when it is perturbed away from a static solu- 
tion. On the other hand, [18] does show that the neglect 



of the back reaction is justified for the parameter values 
which most closely resemble the Goldbcrger-Wise model. 

One might at first feel uneasy about using a 4-D ef- 
fective description of the problem when in reality our 
initial condition is a universe with an infinitely large ex- 
tra dimension. In the Randall-Sundrum scenario, how- 
ever, this is justified because the graviton is trapped on 
the transverse length scale of 1/k <~ 1/M p , rather than 
the size of the 5th dimension (see also [19]). Moreover 
the 4-D Fricdmann equation (63) was shown by ref. [7] 
to be valid without actually assuming the radion to be 
stabilized. § Difficulties with the "wrong" rate of expan- 
sion (H oc p instead of H oc ^fp [20]) arise only when one 
fine-tunes the brane energy densities to prevent radion 
motion even in the absence of stabilization, which we are 
not doing here. In any case, changing the form of the 
expansion rate would have a small effect on our results 
since this alters only the logarithmic term in (65), not 
the overwhelmingly dominant term S. 

The problem of shallow barriers in moduli potentials is 
not unique to the Randall-Sundrum scenario, and a new 
idea for addressing it was recently presented in ref. [21]. 
The coupling of the kinetic terms of matter fields to the 
modulus can give the damping of the modulus motion 
needed to make it settle in the true minimum in some 
cases. This effect might provide an alternative to the 
thermal mechanism we have discussed here. 

We thank Mark Wise and Rob Myers for enlighten- 
ing discussions, Guy Moore for helpful correspondence 
and Genevieve Boisvert for perceptive criticisms of the 
manuscript. 
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